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Abstract 

We study the singular limit of a rotating compressible fluid described by a scaled barotropic 
Navier-Stokes system, where the Rossby number = e, the Mach number = e™, the Reynolds 
number = and the Froude number = are proportional to a small parameter e — t- 0. 
The inviscid planar Euler system is identified as the limit problem. The proof is based on 
the application of the method of relative entropies and careful analysis of oscillatory integrals 
describing the propagation of Rossby-acoustic waves. 

1 Introduction 

The basic system of equations considered in this paper are the continuity equation and the momentum 
equation describing the time evolution of the mass density q = g{t, x) and the (relative) velocity 
u = u{t, x) of a compressible, rotating fluid: 
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dtQ + diY,{gu) = 0, (1.1) 
dt{gu) + div,(pu ® u) + ^^(f X u) + -l-v,p(^) = e"div,s(V,u) + -^gV^G, (1.2) 



§(Va;u) = jj, VxU + V^.u div^ul + ?7div3.ul, > 0, r] >0. 



:i.3) 



The fluid is confined to an infinite slab 

n = R^x (0, 1), (1.4) 

where it satisfies the slip condition 

u ■ n\9Q = [§(V^u) ■ njtanlen = (1.5) 

imposed on the horizontal boundary. 

The model (11.11 - 11.51) may be viewed as a crude approximation (/-plane model) of the Earth 
atmosphere in a plane tangent to the Earth at a certain latitude, see Vallis [Ml Chapter 2, Section 
2.3]. Accordingly, the gravitational force is taken parallel to the vertical projection of the rotation 
axis: 

f= [0,0,1], V,G= [0,0,-1], 

The momentum equation (11. 2p contains a small parameter e related to different characteristic 
numbers resulting from the scale analysis: Rossby number = e, Mach number = e™, Reynolds 
number = Froude number = e^, see Klein [20]. We consider the singular limit problem for 
e \ in the multiscale regime: 

m 

— >n>l, a>0 (1.6) 

for the ill-prepared initial data 

f?(0, ■) = go,e = & + e^g^^K u(0, ■) = uq, (1.7) 

where is a solution to the static problem 

V.p(^,) = £^(™~")^,V,G in n. (1.8) 

In particular, since n > 1, the centrifugal force, parallel to the vertical axis, is dominated by gravi- 
tation (cf. Durran [7]). 
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Formally, it is not difficult to identify the limit problem. Indeed fast rotation is expected to 
eliminate the vertical motion, the vanishing viscosity (high Reynolds number) makes the limit system 
inviscid (hyperbolic), while the low Mach number regime drives the fluid to incompressibility. The 
limit problem is therefore expected to be the incompressible Euler system for the planar velocity 
field V = [vi, V2], 

dtv + div^.(v (g) v) + V^n = 0, div^v = in (0, T) x R'^. (1.9) 

Our main goal is to put these formal arguments on rigorous grounds. 

The phenomena discussed above have been investigated by many authors. The fact that highly 
rotating fluids become planar (two-dimensional), and, accordingly, fast rotation has a regularizing 
effect, was observed by Babin, Mahalov, and Nicolaenko ^J, |2], Bresch, Desjardins, and Gerard- 
Varet, [3], Chemin et al. |4j, among others. The inviscid limit is a well studied and partially still 
open challenging problem, see Clopeau, Mikelic, Robert [5], Kato [I7j, Masmoudi [23], [M], [25] , 
Sammartino and Cafiisch [27], [28], Swann [3l], Temam and Wang [32], [33], to name only a few. 
The low Mach number limits were proposed in the pioneering papers by Ebin [8], and Klainerman 
and Majda [T9], and later reexamined in the context of weak solutions by Lions and Masmoudi [22], 
see also the survey by Danchin [6], Gallagher [H], and Schochet [29]. To the best of our knowledge, 
the simultaneous effects of these three mechanisms has not yet been treated in the literature. 

The present paper may be viewed as complementary to our previous study [Tl], where we exam- 
ined the "single-scale" limit corresponding to 

n = 0, m = l, a>0. (1.10) 

Although this problem looks formally very similar to the present setting, the methods employed as 
well as the limit system are different, cf. [11]. The central issue to be discussed is the behavior of the 
oscillatory part of solutions to the scaled system. These are described in the case f ll.lOp by a system 
of linear equations with constant coefficients, while the situation (11. 6p gives rise to a problem with 
coefficients depending on the scaling parameter e. 

Similarly to [11] , our approach is based on the concept of finite energy weak solutions satisfying the 
relative entropy inequality identified in [13], see Section [2l After collecting the necessary preliminary 
material, we state the main result in Section [31 Section |4] reviews the basic estimates, independent 
of e \ 0, for solutions of the family of scaled problems. Section [5] is the heart of the paper. We 
establish decay estimates on the oscillatory part of solutions to the scaled problems by means of 
careful analysis of certain oscillatory integrals. Here, inspired by the analysis of Guo, Peng, and 
Wang [T6], we make use of frequency cut-off operators and estimates based on van Corput's lemma. 
The proof of convergence towards the limit system is completed in Section [61 
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2 Preliminaries, weak solutions, relative entropy inequality 

We suppose that the pressure p G C[0, oo) fl C^(0, oo) is a given function of the density enjoying the 
following property 

r)'(o) 3 

p{0) = 0, p{g) > for all ^ > 0, Jim 1 > ^- (2.1) 

In addition, without loss of generality, we assume that p is "normalized" setting 

p'{l) = 1. (2.2) 
Finally, we introduce the pressure potential H, 

H{0) = gr^dz, (2.3) 

noting that 

E"{q) = H"(l) = 1. 

e 

2.1 Static solutions 

As already mentioned above, the static solutions solve the problem (11. 8p . specifically, we take 

H'{g,) = e^C™-")^ + H'{1), where G{x) = -Xg; (2.4) 

whence 

geix) = geix3), sup |&(X3)-1| <ce2("-"). (2.5) 

As indicated by our choice of the initial data (11.71) . the solutions of the evolutionary problem (11.11 
- 11.31) . (11.51) . (11.71) satisfy far field conditions in the form 

g ^ Qe, u — 7- as |x| — )■ oo. (2.6) 

2.2 Finite energy weak solutions 

We say that [g, u] is a finite energy weak solution of the problem (11.11 - [L3!) . (11.51) . (11.71) . (12.61) on the 
space-time cylinder (0,T) x i7 if the following holds: 

• Regularity. The functions g, u belong to the class 

> 0, (g- ge) e L~(0, T; + L^n)), u G L''{0, T; W^'^il; R')), u • n = Wglan = 0. (2.7) 



Equations. The equation of continuity f 1 1.1 1) and the momentum equation fll.2p are replaced 
by integral identities 

/ / {Qdt(f + ^'u ■ dx dt = - [ Qo,e'f{0, ■) dx (2.8) 

JO Jn Jn 

for any G C^{[0, T) xU), and 

^ (^QeU^ ■ dt^p + {qu (g) u) : V^v^ - ^g{i x u)-(p + ^p{Q)div^(f^ dx dt (2.9) 

for any G C^{[0, T) xU;R^),^- n|an = 0. 
Energy. The energy inequality 



■) dx + e" [ [ §(V^u) : V^u dx dt 
Jo Jn 

(2.10) 



< 



2^'0,£|uo,e|^ + iH{go^e) - H'{Qe){QQ,e " Qe) " ^(&)) 



dx 



holds for a.a. r G (0,T). 



Note that the existence theory in the class of finite energy weak solutions was developed by Lions 
2T| and later extended in [I2] to the sofar "critical" adiabatic exponent 7 > |- 



2.3 Relative entropy 



For future analysis, it is convenient to replace the energy inequality (12.101) by the relative entropy 
inequality containing more transparent piece of information on the asymptotic behavior of solutions 
for e — 0. To this end, we introduce the relative entropy functional 



r,U 



-^|u - Up + -^(H{q) - H'{T){e - r) - H{t] 



dx, 



(2.11) 



cf. [TO], [13], Germain [15]. 

It can be shown, see [10], that any finite energy weak solution [^), u] specified in Section 
satisfies the relative entropy inequality: 



r,U (T)+e 



" [ In (^(^^") " ^(^-U)) • (V.u - V,.U) dx dt < (2.12) 



^e{Qo,e,Uo,s r(0,-),U(0, 

+ r [ Q {dtV + u ■ V:,.U) • (U - u) dx dt 
Jo Jn 

/ SfV^U) : VJU -u) dxdt + - r f g(i x u) ■ (U - u) da; dt 

JO Jn e Jo Jn 

+4r r / I^"" ~ ^?)^*^'(^) + ^-{H'{r) - H'{~g,)) ■ (rU - gu)] dx dt 
div,u(p(f?) - p(r)) dx dt-^ Jo ""^^-^ ■ U dx dt 



e 

for all sufficiently smooth "test functions" r, U satisfying 

U ■ n\gQ = 0, r > 0, U, (r — compactly supported in Q. 

Note the assumptions concerning the decay and regularity can be relaxed to basically any couple 
[r, U] for which (12.121) makes sense via the standard density argument, see [10]. 

3 Main result 

For a vector field b = [bi,b2,b2], we introduce the horizontal component hh = [&i,&2] writing b = 
[b/ijfes]. Similarly, we use the symbols V^, divh to denote the differential operators acting on the 
horizontal variables only. Finally, the symbol H denotes the Helmholtz projection onto the space of 
solenoidal functions in Q, while H.^ is the Helmholtz projection in R^. 
Let 

vo G W'^^^R^; R^), k>3, divvo = 0, 

be given. It is well known (see for instance Kato and Lai [18j) that the Euler system (II. 9p supple- 
mented with the initial data 

v(0) = Vo 

admits a regular solution (v, 11), unique in the class 

V G C([0, T]; W^'\R^; R^)), dtv G C([0, T]; W''-''^{R^; i?')), U G C([0, T]; W'''^{R^)). 

We are ready to formulate our main result. 
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Theorem 3.1 Let the pressure p = p{g) satisfy the hypotheses ^2. ij) . ( f^.^l) . Suppose that the 
exponents a, m, n are given such that 

m 

a > 0, — > n> I. 

Let the initial data [fpo.e^uo^e] he given by where the stationary states satisfy \2.4\) , 



t'SlU^nL-m) < c, oil -> oi'^ m L^(fi), Uq,,, Uq m R^), (3.1) 



^« e W^-^''^{Sl), uo e for a certain k > 3. (3.2) 

Let [QsyUs] be a finite energy weak solution of the problem ( li.il - UT^) . Iil.5\) . ( [i. 7| j, ( f^.g|) m i/ie 
space-time cylinder (0,T) x f2. 

r/ien 

ess sup \\g^{t, ■) - ^£||(L2+i7)(f^) < e'^c (3.3) 
ie(o,r) 

r weakly-(*) m L'^{0, T; L\n; i?^)), 
v^u, ^ V i (3.4) 
[ strongly m Ll^{{0,T) x n;R^), 

where v = [v/i, 0] is the unique solution of the Euler system ( fi.^) . with the initial data 



V(0,-) = H;, 



The rest of the paper is devoted to the proof of Theorem 13.11 

4 Uniform bounds 

We start with the nowadays standard estimates that follow directly from the energy inequality fl2.10p . 
These are conveniently formulated in terms of the essential and residual components of a measurable 
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function h, 

h = h -\- h 

"' ""ess ~ "res; 

^ess = x{.Qe)h, X ^ C'r°(05 ^)) < X < 1, X = 1 ill ail open interval contaning 1, 

/ires = (1 - 

Since the initial data are given by (11.71) . where the functions Qq\, uo,^ satisfy the hypotheses (13. ip . 
(13. 2p . the integral on the right-hand side of (12. 2p is bounded uniformly for £ \ 0. As the stationary 
states Qe are chosen to satisfy (12. 5p . we deduce the following bounds independent of e\ 



ess sup ||A/&U£||^2m.J^3^ < c, 
te{o,T) ^ ' ' 



ess sup 

te(o,T) 





'Qe - Qe' 








ess 



< c, 



L2(Q) 



ess sup ||[^'e]rcs|ll7(f7) +eSS sup II [1] res 1 1 LI (n) 



and 



iG(0,T) iG(0,T) 

^ /■ 2 

/ VxUe + Va;U* - -div^u^I 

O 



dx dt < c, 



cf. [SI Section 2]. 

Obviously, the estimates (14.11) . (14.21) yield (13. 3p . which, combined with (12.51) gives rise to 

^ 1 in L°°(0, T; L\^J^^)) for any 1 < g < 7. 
Moreover, combining (14. 5 p with (14.11 - 14.30 we obtain 

V^Ue -> u weakly- (*) in L°°(0, T; ^^(fi; R^))^ 

and 

^ u weakly- (*) in L~(0,T; + L^^I^^+^\^- R^)), 

passing to suitable subsequences as the case may be. 

Finally, we may let e — )■ in the equation of continuity to deduce that 



(4.1) 

(4.2) 
(4.3) 

(4.4) 



(4.5) 

(4.6) 
(4.7) 



diVa;U = 0, u ■ n\QQ = in the sense of distributions in (0, T) x Q. 



(4.8) 
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5 Decay estimates and oscillatory integrals 

With our convention (12.21) . the equation describing the oscillatory part of solutions reads 

e'^dts + div^V = 0, (5.1) 

e'^dtV + wf X V + V^s = 0, 00 = e'^-\ V ■ n\aa = 0, (5.2) 
cf. [11]. Re-scaling in the time we arrive at 

dts + div^V = 0, (5.3) 

dtV + cuhxV + V^s = 0, V ■ n\9n = 0, (5.4) 

with the operator 

Kf \ \ ^ ^ .\ div^.V 
Biuj : t-). 
^ ^ V cjf X V + 

The operator B is skew symmetric in the space L'^{Q) x L^(fi; R^), with the domain of definition 
V[B{oo)] = {[r, V] I r G W''\n), V G L^Q; R'), div.V G ^^(n), V ■ n = ^gbn = o} , 
and the kernel 

Af{B{u)) = {[g, v] I g = g(x;,), q G iy^'2(i?2), v = [v;,(a;;,), 0], div^^v;, = 0, x v + V,g = o} . 
Let V{ijj) denote the projection 

V{u) : L'^iQ) X L\Q;R^) ^ Ar{B{u)). 
Exactly as in [lU Section 4.1.1] we can show that 

P(a;)[r,U] = [g,v] 

if 

-Ahq+u'^q = u I cmlhVh dxs+w^ / r dxs in i?^ v = [vi, V2, 0], uvi = -d^^q, UV2 = d^^q. (5.5) 
JO Jo 
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5.1 Spectral analysis 

Thanks to our special choice of the geometry of the spatial domain Q, we may reformulate the 
problem (15. 3p . (15. 4p in terms of the Fourier variables. To this end, we observe, exactly as in 
that the underlying spatial domain Q may be equivalently replaced by 

where 

r^ = [-i,i]{-i,i} 

is a "flat" sphere, and where s, V/j were extended as even functions of the vertical variable 0:3, while 
V3 was extended as odd in x^. In other words, all quantities are understood as 2— periodic functions 
in the vertical X3 variable. 

Accordingly, for each function g G L^(fi), we introduce its Fourier representation 

di^y = —?K I / exp (-i^ ■ Xh) g{xh, X3) exp (-i/cxg) dxh dajg. 
V2 J-i Jr2 



where 



We have 



g{xh, X3) = ^^Xxn [9{i, k)] exp (i/cxg) 



where the symbol J^x^^i denotes the standard Fourier transform on B?. 
Thus the problem (15.31) . (15.41) can be written in the form 









d 




+ i 


It 






V^{t,i,k) 





" 


^1 


6 


k ' 











(jji 







Vi{t,^,k) 


6 


—ui 










V2{t,^,k) 


k 








. 




Vsity^yk) 







So{^,k) 






. Vo(e,fc). 



(5.6) 



with the Hermitian matrix 



A{^,k,u;) 






6 


6 


k 







ui 







—ui 








k 












Of course, solutions of (15.61) depend also on the parameter to = e 



m—l 
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5.1.1 Spectral properties of the matrix A 

After a bit tedious but straightforward manipulation (see [9]), we can check that A{^, k,u) possesses 
four eigenvalues 



w2 + |^|2+fc2+^(a;2+|5|2+fc2)2-4a;2fc2 
2 



i^2 + |5|2+fc2_^(^2 + |g|2+fc2)2_4^2fc2 



1/2 



1/2 



A2(|eP,fc,a;) = -Ai(|ep,fc,a;), 



A4(|eP,A^,c^) = -A3(|eP,fc,a;). 



(5.7) 



Note that A3(|^p,0,a;) = A4(|^p,0,a;) = are the zero eigenvalues corresponding to the non-void 
kernel of the matrix A- the Fourier image of the null-space of the operator B{uj), see [9]. 
As for the eigenvectors [q, f i, f2, fs], we have 

^iVi + ^2V2 + kvs = Xq, ^iq + iuv2 = Xvi, ^2q - = Xv2, kq = Xvs, (5.8) 

from which we immediately deduce 

vi = iJ,{X^i + iw^a), V2 = /i(A^2 - iw^i), /iA|^p = Ag - kv^, , kq = Xv^, (5.9) 

where /x is a free parameter that is fixed to normalize the norm of the eigenvector to be one. 



5.1.2 Eigenvectors 

We denote by E = Fi{^,k,ijj) = [q,Vi,V2,V3,] the normalized eigenvectors. Our goal is to show 
that diagonalizing matrices Q, Q^, formed by the eigenvectors, are L^— multipliers in the ^ variable 
restricted to compact subsets of \ {0}. This amounts to showing that 



sup max 

j6(0,1) 0<a<|^|<6<oo 



VfEji^, k,uj) <c = c{A, a, 6, k), j = 1, 2, 3, 4, ke Z 



(5.10) 



with any multi- index A = {Ai, A2). 
We distinguish two cases: 

Case A = 0: 

In this case, we necessarily have k = and A = iAs, and the orthonormal basis of eigenvectors 
can be taken in the form 



El = /i 
Clearly flCTl holds. 



6,6,0], /x= + 



co'') E2 = [0,0,0,1]. 
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Case A ^ : 

We find that 



A 



Thus, the corresponding normahzed eigenvector has the form 

E = /i 

with 



A^i + A6 - iw^i, -^^^ 

-1/2 



A4+ |A;|2A2 



.ier+(A^+u;^)ier 



(A2 - F)2 

We consider first the case A = ±Ai. We check that 

A?>ier/2, {u'+\e+kr-4u'e>\^\\ 



{5.11] 



Consequently flS.lOp is satisfied. 

Finally, if 7^ A = iAs, we note that 

x^m\k,u)=u\k\ 



Ai(iep,fc,^: 



Xl-k^ 



k\uj^ - \l) 



' ■• — > 1^1 /2. 



(5.12) 



Using identity fl5.12l) i . we write E in terms of Ai and verify (15.101) employing this explicit formula 
and estimates (IKTTTll . f l57T2l) . 

5.2 Frequency cut-off 



As we shall see below, it will be convenient to approximate the initial data for the problem (15. 3p . 
(15. 4p by a frequency truncation represented by a function ip G C^(0, 00). 

Accordingly, solutions of the problem (15. 3p . (15.40 will be, in a certain way, composed of the 
quantities 

Z{t,Xh,k,u) = J'^X., [exp(±iA,(|e|2,A;,a;)t)^(|e|)/i(0] , J = 1, 2, 3, 4, e Z, G (0,1), (5.13) 
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where h stands for the Fourier transform of the "initial data". Our goal will be to derive suitable 
dispersive estimates for the mapping h Z. 

We start with the — L°° decay estimates. To this end, write 



\\Z{t,-,k,Uj)\\Loo(^R2) < 

Consequently, 



\h\\Ti 



TT 



exp(±iA,(|eP,fc,u;)t)^(|e|y 

J-^-^,^ [exp(±iA,(|er,fc,u;)t)V^(|el)] K) 
exp ( ± iA,(|eP, uj)t)i^m) exp (i|elk.| sin^) \i\ ^dO 



(5.14) 
(5.15) 



Jo Jo 



7r 



V2 



exp 



( ± iAj(r^, k, uj)t^ilj{r)rJo{r\xh\) dr. 



where the symbols Jm, m = 0, 1, . . . denote the Bessel functions, cf. Guo, Peng, and Wang [16]. 
Finally, performing a simple change of variables, we get 



exp ±iA,(|ehfc,a;)th^(|e|) 



K^h) = -^r- I exp ( ± iXj{z, k, u)t)'^p{^/z)Jo{^/z\xh\) dz. 



2 Jo 



(5.16) 



5.3 Decay estimates 



Supposing Xj 7^ we derive the desired decay estimates. To this end, we use van Corput's lemma, 
see Stein [30], Chapter 8.1.2, Proposition 2 and Corollary]: 



Lemma 5.1 Let A = A{z) be a smooth function away from the origin, 

dzA{z) monotone, \dzA{z)\ > Aq > 
for all z E [a,h\, < a < b < oo. Let ^ be a smooth function on [a, b] . 

Then 



exp {iA{z)t) ^(z) dz 



< c 



1 

tTo 



J a 



where c is an absolute constant independent of A and $. 
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Going back to the oscillatory integral fl5.16p . we distinguish two cases. 



5.3.1 Case \xh\ > 

Using the decay properties of Jo and the fact that ip is compactly supported away from zero, we get 
exp(±iAj-(|^P,A;,w)t)7/;(|^|)] < c(V^)r^/^ whenever \xh\ > ■ (5.17) 



5.3.2 Case \xh\ < 

The idea is to use van Corput's lemma. Let [a, h] be a closed interval, a > 0, containing the support 
of iply/z). In accordance with the hypotheses of Lemma [5. H we have to verify that 



\dzXj{z, k,ijj)\ > Ao(a, b, oj,k) > for all z G [a, b], uj G (0, 1); 



• <9^Aj is monotone on [a, 6]. 

If the two conditions are satisfied, we get 



exp 



1 



(l + t'^/^^ for Ix;,! <t^, (5.18) 



Ao(a, b, k, uj)t 

where again we have used the properties of the Bessel functions, namely, 

J'oiz) = -Ji{z). 

Now, our goal is to verify the hypotheses of Lemma 15.11 We have 

d;,Xi{z,k,u) = ^ , , -(3i{z,k,uj), (3i{z,k,uj) = ^ ( 1 + ^ +z + k 

2\i{z,k,uj) 2y J(u;2 + z + F)2_4a;2^2^ 



and, similarly. 



1 1 



2 Xs{z, k,u) 



dAsiz, k, oj) = ^TT— 7— t/^3(^, k, u), Psiz, k,u) = ^\ l " ' ^" ) , k ^ 0. 



u'^ + z + k"^ 



'uj^ + z + k^f-Auj^k^ 



Furthermore, 

d^l3i{z,k,u) 



2uj'k 



2u2 



[(u;2 + ^ + A;2)2-4a;2^2]3/2' 

14 



d^(5^{z, k,u) 



2u2 



[(^2 + P)2 - 4a;2^ 



21,213/2 ■ 



Summing up the previous relations, we conclude that 

dzXi{z, k,uj) is a decreasing functions of z, 

while 

dzXsiz, k, u) is an increasing functions of z for k ^ 0. 
Consequently, we deduce that 

d,Xi{z,k,u) > d,Ai{b,k,u) > A{^,k) for z G [a,b], u e (0,1). (5.19) 

Finally, we have 

dzX3{z,k,u)) < 0, \dA3{z,k,u)\ > {dAsib, k,u)\, 

where 

-1/2 



I co' + b + e \ fu' + b + k'-^{u^ + b + k^y-4co^k^ 

\d,X3{b,k,uj)\ > - - 1 



= \/2uj\k\ ((a;2 + 5 + _ 4^2^2^ 1/2 |^^2 ^ ^ ^ ^2 ^ y^(a;2 + 5 + p)2 _ 4^2^ 

> c{^,k)uj for A; ^ 0, we (0,1). 
Thus, reviewing f l5.17p . (15.181) we may infer that 



-1/2 



•^^4.. [exp ( ± iA,(|eP, fc, a;)t)^(|e|)] ||^^(^,^ < c(V^, fc) max {-^^; ^} , t > 0, (5.20) 
as soon as Xj 7^ 0, which gives rise to the decay estimates 

\\Z{t,-,k,u)\\L^(^R2) < c{ij,k)maxl^—^^; ^| \\h\\Li(R2). (5.21) 

Next, seeing that the mapping h exp {iXj{^, k,u)tjh is an isometry on L^(i?^), we deduce 

\\Z{t,;k,u)\\L2iR2)<c\\h\\mR2) (5.22) 

Finally, interpolating f l5.2ip and (15.221) . we obtain the — L'^ estimates 

f 1 1 1 ^"f 11 

\\Z{t,-,k,u)\\LP(R2) < c(7/;,p,fc)max| ^^^_^^^ ;^| for p > 2, - + - = 1, ^ 0. 



p p 

(5.23) 
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Keeping in mind that u scales like ^ while the time t is proportional to e we observe that 
taking 

2 

< /3 < — 
m 

yields the effective decay of = Z{t/e"^,k,u) on any compact subinterval of (0,T]. In particular, 
the optimal choice /3 = 1/m gives rise to 



^ i_i I 1 1 I^-P 

— ^1— l/2m ' ^l/2m 



/^ILp'(fi2), p>2, 7^0, t G (0,T]. 



(5.24) 



6 Convergence 



In this final part, we complete the proof of Theorem 13. II The basic idea is to use the relative entropy 
inequality (12.121) for a suitable choice of test functions r and U. 



6.1 Initial data decomposition 

We start be introducing suitable smoothing operators imposed on the initial data. Taking a family 
of smooth functions 

ijjs e C^{0, oo), < ^5 < 1, / 1 as (5 ^ 0, 

and 

<Ps = (psM e C^{R^), < 05 < 1, 05 / 1 as 5 ^ 0, 

we introduce 

1 



exp {—ikx3 



So,e,5 + qo,e,5, where - Ahqo,e,5 + (^^qo,s = w / curlft [[no]h]s dxg + u 



1 r 



(6.1) 



and, similarly, 

\^^,]\A^h,^z) = ^ E [V^5(l^l)(wo,j05)(^,/i;)] exp(-iA;x3) , j = 1,2,3. (6.2) 

Now, we write the initial data in the form 



.dx3 (6.3) 



[Uo]^ = Vo,£,5 + Vo,£,5, with Uj\vQ^eA^ = -9x290,8,5, Uj[vo,e,s]2 = dx^qo^e,5- 



(6.4) 
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Finally, we choose the functions r, U in the relative entropy inequality as 

f = re,5 = 06+ £"'iqe,5 + Se,^), U = \Je,5 = Ve,5 + ^e,5, (6.5) 

where [s^^^, V^ ^] is the unique solution of the acoustic- Rossby system (15. ip . (15.21) . emanating from 
the initial data 

5^,5(0, ■) = So,e,5, Vs,5(0, ■) = Vo,e,S, 

while the functions q^^s, ^e,5 are interrelated through 

oof X Vs,5 + V^qe,5 = 0, (6.6) 

where qs^s is the unique solution of the problem 



dt (^hqe,5 - UJ^qs,5) + ^"^^(16,6 ■ (^Ahqe,5 " ^^^e.^) = 0, ^^,5(0, ■) 



(6.7) 



6.2 Decay of the oscillatory component 

First we claim that, in view of the dispersive estimates (I5.24p (with < /3 < 2/m), and the properties 
of the eigenvectors of the matrix A, discussed in detail in Section I5.1.2[ we get 



Se,s ^ in 1^(0, T; W^'\n)), Vs,5 ^ in 1^(0, T; W^^'^{n)) as £ ^ 



(6.^ 



for any fixed 5>0, l<p<oo, 2<g<oo and Z = 0, 1, 
We note that 



S£,s(t, Xh, k) 



exp(iAi(e,a;,A;)^),0,0,0 

0,exp (iA2(^,a;,/c)^) ,0,0 

0,0, exp (iA3(^,a;,/c)^) ,0 

0,0,0,exp {\\i{^,uj,k)^ 
exp (iAi(^,w, A;)^) ,0,0,0 



Q(e,a;,A;)^(|e|)/^o(0 



0,exp (i\2{i,u,k)^^ ,0,0 
0,0, exp {i\^{i,uj,k)^) ,0 
_ 0,0,0, exp {i\i{^,uj,k)- 
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m\)Q{i.^.k)m\)ho{i) 



u = e 



.m—l 



where ip G C^(0, oo) has been chosen so that ipip = 

In accordance with Section [5.1.21 the quantities tpQ, tpQ'^ are L^— Fourier multiphers with norm 
independent of u. Thus the desired decay estimates (16.81) follow from (15.241) . 

6.3 Convergence of the non-oscillatory component 

We introduce the scaled function (jir^s = <le,s/<-^ and observe that 

dt [Ahqs,s - oj'^qe,&) + ^iqe,s ■ Vft [Ahqe,s - ^'^qe,s) , 9e,5(0) = go,M' (6-9) 

Ahqo,e,5 - (^^qo,5 = f curl/, [[Uo] ft] 5 dXg + / LJ,^^] dXg. 

Jo JO L J5 

Since the initial data are regular, we may use the result of Oliver f2E[ Theorem 3] to deduce that 

{Aq.^s + £^qe,5}e>o IS bouuded in C''([0,T]; W^^^{n)), r > 0, / = 0, 1, . . . (6.10) 

Moreover, as Aq^^s satisfies that transport equation (16.91) with the initial data in L^(i?^), we get 

{Aq.^s + uJ^qe,5}e>o IS bouudcd in L°°([0, T]; LP{Q)) for any 1 < j9 < cx). (6.11) 

Next, we recall the "energy estimates" that can be obtained multiplying (16.91) on q^^s and inte- 
grating by parts: 

J^^ {\^xqe,s\^ + oo^\qe,s\^) (r, ■) dx = J^^ {\V^qo,e,s\^ + oo^\qo,e,s\^) dx. (6.12) 

Note that 

/ Vft ■ VhAhqe,5qe,sdx = - ^iqe.s ' V hqe,5Ahqe.sdx = 0. 



Since 
we get 

K,5}e>o is bounded in C"^([0,T];iy''2(i?2))^ ^ > 0, / = 0,1,... (6.13) 
Finally, we compute 

dtqe,s = {Ah - u^y^ [Vh {ve,scmWe,s)] , 

where, furthermore, 

vcurlftV = vi {d^^V2 - d^^vi) ; V2 {d^iV2 - d^^Vi) 



with 
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Consequently, we may infer that 

{dtqe,s}e>o is bounded in C^([0,T]; W^'^iR^)), r > 0, g > 1, / = 0, 1, . . . (6.14) 

All the above estimates may depend on 6 but are uniform with respect to e \ 0. 

In view of the above estimates, it is easy to pass to the limit for e — )■ in order to get 

^e,5 ^ V5, 9iV,,5 ^ ^t^rs weakly-(*) in L°°(0, T; W'^^R^)), 1 = 0,1,..., (6.15) 

where 

V,, dtVseCi[0,T];W''\R^)), / = 0,1,2,.... (6.16) 

We have, in particular, 

v,,5 -> V5 in L%0, T; Lf„,(fi)), 1 < q < oo; 

whence, by virtue of (16. 9p . 

dtcmlh^s + V5 ■ VhCurUv^ = 0. 

Seeing that 

curl/,div/j(v5 (g) v^) = curU(v5 ■ V^v^) = ■ V^curUv^, 
we deduce existence of 

nseC{[0,T];W'^\R'')), / = 0,1,2,..., (6.17) 

where the couple (v^. Us) is the unique solution of the Euler system (II. 9p . emanating from the initial 
data 

V5(0,-)= Th^H^] dx3. (6.18) 





6.4 Relative entropy inequality 

We return to the relative entropy inequality (12. 12^ . where g = Qs, u = and the test functions r 
and U are given by (16.51) . In what follows, we examine step by step all terms on the right-hand side 
of (I2.12P and perform the limits; first for e — )• 0, then for 5 — )■ 0. 
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6.4.1 Initial data 

We have 



r(0,-),U(0, 



(6.19) 



Jn I 



< c(||uo,e - [uo]5||i2(f^.^3) + Wq^^I - [^0^^]5||i2(o.^3)) 

^ c (||uo - [yiM\h(n;R^) + 11^0^^ - [S\\\h{n;R:^)) as e ^ 0. 
The most left quantity obviously tends to zero as 5 \ 0. 

6.4.2 Viscosity 

We write 



(i)i 



10 Jn 

and, by Korn's inequality, 



S(V,U,,5) : V,(U,,5 - u,) dx dt 



<£"ci((5) / / |V,(U,,5-u,)| da;dt, 
JO Jn 



(6.20) 



e" / / |V:,(U,,5 - Ue)| dx dt < — / / (Se(V^Ue) - Se(V„Ue,5) : (U,,5 - u,) dx dt + cs^". 
/o Jn 1 Jo Jn 



6.4.3 Forcing term 

We have 



1 



r2m 



(6.21) 



+ ^div^Ue,5(p(r^,5) -p(&)) + -^{r.^s - Qe)^xG ■ lJe,5 



r2m 



pir^^s) - p'{re,s){re,6 - &) - p(&) div^Ue,^ + -^{re,s - ge)H"{re,s) dtr^^s + div^{re,s^e,6 



r2m 



QeV^H\re,s) ■ (U,,5 - U,)-V.iy'(&) ■ {n^S^e^S ' feU,)] + -J^{re,5 ' &)V,.G • U,, 



where 



dtr^^s + d\v^{r^^s^^^s) = (& - l)div^V£_5 + V^^^^^ ■ 11^,5 + 

e"'dtqe,5 + e'"div^. ((Se,5 + qe,5){^e,5 + V,,^)) • 



(6.22) 
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and 

p,V,i/'(r,,5) ■ (Ue,5 - u,) - V^H'ig,) ■ (r.^^U,,^ - ^,u,) (6.23) 
= gsV.[H'{r,,s) - H"{~g,){r,,s - Qs) - H'{g,)] ■ (U,,^ - u,) 
+0s^xH"{ge) ■ (U,,5 - u,)(r,,5 - g,) + gee"'H"{g,)(y^s,,s + ^ .qe^s) ■ {^e,5 - u,) 

HQs - re,s)V,H'{g,) ■ U,,^. 

Here we have used the fact that the term ge^ xH'^g^) ■ (U^^^ — u^) in the expansion of ge^ xH'{r^^s) ■ 
{\Je,5 ~ Ue) cancels with the same term in the expansion of VxH'{g^) ■ {r^^s^^^s — Qe'^e)- 

Now we use formula fl6.2ip with the second and third terms at the right hand side expressed 
through formulas (I6.22p and (16.231) : we keep the terms ^gsH"{gs){yxS£,5 + '\^xC[e,5) " (U^^^ — Ug) and 
^(^^£,.5 - Qe)^xH'{r^^s) [dtqe,5 + diVa; ((5^,5 + q£,5){^e,5 + ^e,&))] as they are, and estimate the decay of 
all remaining terms as e — )■ 0. 

We observe that, thanks to ([63]), fl6:T2|l . 



sup \\qe,5it,-)\\L^R^) <c, sup || Ag^,5(t, ■) ||i2(^2) < ce™ , sup \\dtqe,sit,-)\\L^(R^) < ce 
t€{(},T) te(o,T) ie(o,r) 



m—l 



Now, we use (EaHZSl), fOM . fl63D . (KSh . (Km . (EMD to deduce the following estimates: 



(6.24) 



1 



r I {re,S - Qe)VxG ■ U,,5 dxdt 



< ce 



m—2n 



(6.25) 



£2™ Jo Jn 



e^™ Jo Jn 



P{rs,5) - p'ire,5){re,5 - Qe) " P(&) 
pir-e^s) -p{re,s){re,S - ge) " p(&) 



divT.Uf A dxdt 



div^Ve ,5 dxdt 



£2"^ Jo Jn 



Pire,5) - p'ir£^s){rs,5 - Qe) - p{Qe) ciiv^V,,,^ dxdt = h{£, 5) 



(6.26) 



r2m 



Jn 



re,s - Qe)H"{r^^s)({Q6 - l)div^.V£,5 + V^ge ■ Ve,5j dxdt 



< ce 



m—2n 



Jo Jn 



f?,V. H'in^s) - H"{g,){r,,s - Qe) - H\~g,) ■ (U,,^ - u,) dxdt = h{e, 5) 



(6.27) 
(6.28) 
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^( f geV^H"ig,) ■ (U,,5 - u,)(r,,5 - g,) dxdt + T f - r,,5)V^H'{g,) ■ U,,^ dxdt 

Here and hereafter, (5) denotes a generic function having the property 

h{e, 6) -> h{6) a.s e ^ 0, h{5) as 5 0. 



< ce 



m—2n 



(6.29) 
(6.30) 



Note that the dispersive decay estimates (16. 8 p play a crucial role in the analysis. 

Taking into account ( ]6.19p . (I6.20p . using the identity (I6.2ip with the third and fourth terms 
expressed through (I6.22f[6l23l) . and employing the asymptotic behavior from formulas (I6.25f[6l29|) . 
we may rewrite the relative entropy inequality (12.121) in the form 



1 

+ — 



+ - 



< r f ge {dtV.^s + ■ V,U,,5) ■ (U,,5 - U,) dx dt 

Jo Jn 

+ - f I &(f X ■ i\Je,5 - Ue) dx dt 

e JO Jn 

Te^S - ge)H"{rs,5) dtqe,5 + div^ {{Se,5 + ge,^)!!^,^) 

^ &^"(&)(v,s,,5 + V.g,,^) ■ (U,,5 - u,) dx dt + h{e, 5). 



(6.31) 



dx dt 



6.4.4 Coriolis force 

We may write 

geH"{~g,) (V,s,,5 + ^.qe,s) ■ (U,,^ - u,) (6.32) 

= ge{H"{~ge) - H"il))(y,Se,5 + V,q,,s) ■ iUe,5 ' U,) + &/7"(l) (V.S,,5 + Vccqe,s) " (U,,5 - U,), 

where, by the same reasoning as in estimates (I6.25f[09i) . 
1 



jje{H'\~g,)-H"{l)){y,s,^s + V^q,,s) ■ (U,,^ - u,) dxdt 
Recalling our convention H"{1) = 1 we get for the second term in (16.321) . 

■^geH"{l)(y,s,^s + V^g,,^) ■ (U,,^ - u,) = 

- dt\e,5 - -fx V,^s - -fx V^,^) ■ (U.^s - Ue) 



< ce 



m—2n 
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where the last term cancels out with the Coriolis force term -f x ■ (Ue^s — u^) = -(f x u^) ■ Ue^s- 
Thanks to formula f l6.32p the relative entropy inequality f l6.3ip reduces to 



e'"- Jo Jn 



Ue,5) (r) < r f Qe {dtY,,s + u, ■ V.U,,^) ■ (U,,5 - u,) dx dt (6.33) 
^ Jo Jn 

ire,5 - Qe)H"{r^^s) dtqe,5 + div^ ((5^,5 + qe,s)'U^^s) dx dt + h{e, 5). 



dt 



6.4.5 Estimating the remaining terms 

We observe that 

/ / & {dt'^e,S + ■ V^Ue,^) ■ (Ue,5 - U^) dx dt = [ j Qe {^e,5 + U^^^ ■ V^\Je,5) ' {^e,5 " ^e) dx 

Jo Jn Jo Jn 

- [ f ge{'Ue,5 - Ue) ■ VxVe,5 " {'Ue,5 - ^e) dx dt, 

Jo Jn 

where the terms in the first expression at the right hand side containing the quantities s^^s, ^e,s tend 
to in the limit lim^-j^o lim£_>.o thanks to dispersive estimates (16. Sp . Consequently, employing (I2.5p . 
we obtain ^ 

Qe {dtVe.S + Us ■ V^U^.^) ■ (lJe,S - U^) dx dt 

Jo Jn 

~ lo X "^^'"^ ' ' '^^^'■^ ~ dxdt + cj^ Se[ge,Ue rs,s,'Ue,5) dt + h{e,5). 

Similarly, we deduce 

dx dt 



Jo Jn 



{re,5 - Qe)H"{re,5) dtqe,5 + div^ {{Se,5 + qe,5)Ve,5 

= [ [ qe,sdtqe.s dx dt + h{e, 6), 
Jo Jn 

where we have used (ICTn . fle^]) . fICTD . 

Now, we employ the energy equality (16.121) to observe that 

J (ve,5 ■ dtVe,S + qs,e " ^t^e.^) dx dt = 0. 



Consequently, 



/ / {dtVe,5 + Ve,5 ■ V:,Ve,s) " {^e,5 - U^) dxdt 

Jo Jn 

+ f f qe,5dtqe,5 dx dt = - [ [ {dtVe,5 + ^6,5 " V^v^^^) ■ dxdt, 
Jo Jn Jo Jn 
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dt 



where we have used the identity Jq ,5 ■ V^v^ ,5 ■ ,5 dxdt = 0. 
We remark that 

/ / {dtVe,5 + Ve,5 ■ VxVe,s) ■ dx dt ^ [ [ {Ot^s + vs ■ Vx^s) " u dx 
Jo Jn Jo Jn 

= f f VxTls ■ u dx dt = 0, 
JO Jn 

where we have used (I4.7H4.8I) . f l6.15p and equations (11.91) with (v, 11) replaced by (v^,!!^). Conse- 
quently, the entropy inequality f l6.33l) takes the form 

r,^s, ^6,5) (r) < h{e, 6), 
where the function h satisfies (16.301) . This finishes the proof of Theorem 13. II 
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